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The velocity-dependent ‘one-scale’ model of Martins & Shellard is used to study the evolution 
of a cosmic string network (and the corresponding loop population) in open universes. It is shown 
that in this case there is no linear scaling regime and that even though curvature still dominates 
the dynamics, at late times strings become the main component of the universe. We also comment 
on the possible consequences of these results. 
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I. INTRODUCTION 

Despite the strong theoretical prejudices favouring a flat universe [QQ, there is a fair amount of observational data 
which suggests the possibility of an open universe, with a present density that could be as low as O ~ 0.3—notably, 
the so-called ‘age problem’ Q and recent measurements of the baryonic content of x-ray clusters [Q. Since in a 
low-0 universe structures collapse earlier, observations of galaxies at high redshift would also be easier to explain if 
the universe is open. It is therefore appropriate to consider how some of the standard cosmological scenarios would 
change if this possibility turns out to be true. 

One relevant case is that of the evolution of a network of cosmic strings ||. It has been shown HQ that defect 
models normalized to COBE in an open universe predict a galaxy power spectrum consistent with that inferred from 
galaxy surveys without requiring an extreme bias (in general, fl = 1 models predict more small-scale power than 
low -Clh ones). However, these results were established either using a priori scaling assumptions for the string network 
H or numerical simulations of texture evolution [Q. Here we study the evolution of a cosmic string network in 
open universes, using the velocity-dependent ‘one-scale’ model of Martins & Shellard HQ], which provides the first 
quantitative description of the complete evolution of the large-scale properties of a cosmic string network. This is 
briefly summarised below, and used in the following section to obtain the evolutionary properties of both the long- 
string and the loop populations in an open universe; these are then compared with the standard flat universe case, 
and some implications of these results are discussed. 


II. THE VELOCITY-DEPENDENT ‘ONE-SCALE’ MODEL 

Two different but complementary approaches have been used to study cosmic string evolution. The simplest 
(although more expensive) is to use large numerical simulations |r]]. Among other interesting things, these revealed 
a significant amount of small-scale structure (or ‘wiggles’) on the strings, containing up to one half of the total string 
energy. 

On the other hand, there is always thew possibility of using analytic methods—an approach first used by Kibble 
p[ . Due to the strings’ statistical nature, what one really does is ‘string thermodynamics’, that is describing the 
network by a small number of macroscopic (or ‘averaged’) quantities whose evolution equations are derived from the 
microscopic string equations of motion, and introducing additional ‘phenomenological parameters’ if necessary. The 
first such model providing a quantitative picture of the complete evolution of a string network (and the corresponding 
loop population) has been recently developed by Martins & Shellard ||,Q; this has the added advantage of being 
equally applicable to the study of vortex-string evolution in a condensed matter context ||l2|j . Here we will present a 
very brief description of this model—the reader is referred to the original paper [Q for further details. 

Apart from the straightforward definition of the energy of a piece of string, E = iia(r) J eda (e being the coordinate 
energy per unit er), the only other macroscopic quantity in this model is the string RMS velocity, defined by 
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Explicitly distinguishing between long (or ‘infinite’) strings and loops, we can use the fact that the former should 
be Brownian to define the long-string correlation length as p^ = p/L 2 . Phenomenological terms must be included 
for the interchange of energy between long strings and loops. A ‘loop chopping efficiency’ parameter (expected to be 
slightly smaller than unity) is introduced to characterise loop production 
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\dt~ 


to loops 
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and in the particular case of GUT-scale strings (but not more generally |9|) it can be safely assumed that loop recon¬ 
nections onto the long-string network are negligible. In particular, this has been confirmed in numerical simulations 
0- Note that it is conceivable that the behaviour of c is different in flat and open universes. However, this effect 
will not be crucial, because the scaling properties do not depend strongly on it 

It is then simple to derive the evolution equation for the correlation length L. Since we are only interested in the 
epoch around radiation-matter equality, we need not be considering frictional forces Q , and we simply have 


2— =2HL(1+v 2 00 )+cv 00 
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For the case of string loops, the relevant lengthscale is simply the loop length, which decays due to gravitational 
radiation, and its evolution equation is 


df 

- = (l-2v 2 )m-r'Gpvf, (2.4) 

where T' ~ 8 x 65. Then the only other thing that is needed is an assumption on the loop size at formation. In 
the epoch relevant to this paper, we expect it to be approximately constant and much smaller than the correlation 
length—we will take ii = 10 ~ 3 L(ti). Then for any given time, one only has to look at the loops that have formed 
until then, determine which of them are still around and add up their lengths to determine the total energy density 
in the form of loops. This is conveniently expressible in terms of the ratio of the energy densities in loops and long 
strings 


6{t) = —jjr = gcL 2 (t) 
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where g is a Lorentz factor accounting for the fact that loops are usually produced with a non-zero centre-of-mass 
velocity, a is the loop size at formation relative to the correlation length at that time and £(t, t') is the length at time 
t of a loop that was formed at a time t'. 

Finally, one can derive an evolution equation for the long string or loop velocity with only a little more than 
Newton’s second law 


!=( i-” 2 ) 

here k is another phenomenological parameter that is related to the presence of small-scale structure on the strings; 
an appropriate ansatz for it is (refer to || for a complete justification) 

k-{\ 2HR>X (27) 

l 72 2HR ’ 2HR <X 1 l) 

where R is the curvature radius of the string (that is , R = L for long strings, but l = 2 ttR. for loops) and x is a 
numerically determined coefficient of order unity, whose precise value depends on whether one is using the above 
ansatz for long strings or loops—see || for a complete discussion of this point. 

The above quantities are sufficient to quantitatively describe the large-scale characteristics of a cosmic string network 
around the epoch of equal matter and radiation densities (see Q). In a more general situation one would need to 
include the effect of frictional forces |jj due to particle scattering on the strings. 


2 






III. EVOLUTION IN AN OPEN UNIVERSE 


We now come to the issue of this paper. As we already pointed out, we only need to study the behaviour of the 
string network in the transition between the radiation- and the matter-dominated regimes (see Martins & Shellard 
H for a detailed discussion of the early stages of evolution of GUT-scale and other cosmic string networks). 

It is straightforward to see that there is one crucial difference with respect to the case of a flat universe: in an open 
universe there will no longer be a linear scaling regime. This arises naturally from the fact that in a universe where 
the scale factor grows as a oc t x (with A < 1) the linear regime has the following properties 
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In an open universe the ‘effective’ A is a variable, increasing from A = 1/2 in the radiation era to an asymptotic value 
of A = 1 (see figure [j]). For example, if we happen to live in a universe where fl 0 ~ 0.3, we have A 0 ~ 0.8 today, 
and A 0 ~ 0.9 for fl 0 ~ 0.1. In other words, there will be corrections to the simple linear behaviour, such that the 
correlation length L will grow slightly faster than t. Nevertheless, since the horizon size for an a oc t x universe is 

dn = ’ ( 3 - 3 ) 


one can easily show that L will always be smaller than the horizon. On the other hand, the string velocity will 
decrease with time, and therefore loop formation will gradually switch off. 

Note that the power-law dependence of the scale factor obviously changes in the transition between radiation and 
matter domination, even in the case of a flat universe. In particular, there will of course be a departure from linear 
scaling while the transition is taking place, with L/t growing from 0.27 to 0.6 (approximately). 

In particular, we can easily find the solution to the averaged equations of motion in the limit where A = 1, 

L oc t (Ini) 1 / 2 , v oc (Ini) -1 / 2 ; (3-4) 

note that in this limit the correlation length grows more slowly than the horizon (which goes like dn oc tint). 

In order to get quantitative results, we must solve the averaged evolution equations described in Section || numeri¬ 
cally. To these we must add a further equation—the Friedmann equation—specifying how the scale factor (and hence 
the Hubble parameter) evolves in the transition between the radiation- and the matter-dominated epochs: 
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where K = 0 for a flat universe and K = — 1 for an open universe; note that we should at least consider the possibility 
of the string density becoming a non-negligible source for the Friedmann equation. 

Figures ^fand|^ contrast the evolution of the long-string and loop populations in a flat universe and in open universes 
with = 0.3 and S! 0 = 0.1 (we have assumed that h = 0.6 in both cases); note that the present epoch corresponds 
to a 0 /a teq ~ 2.3 x 10 4 Q 0 h 2 . 

As was first discussed in ||, even in the case of a flat universe the transition from the radiation to the matter epoch 
is a very slow process, lasting about eight orders of magnitude in time. In the case of an open universe, apart from the 
differences we already expected, the most interesting result is that, although the string density always decreases with 
respect to the critical density (as one would also expect), at a redshift around z ~ fl” 1 (which is approximately when 
curvature has started to dominate the dynamics) the string density has started to grow relative to the background 
density. In fact, in a fl 0 = 0.3 universe strings will become the main component of the universe in about seven orders 
of magnitude in time, whereas if we had f l 0 = 0.1 this would only take about four orders of magnitude in time. 

It is also interesting to point out that, although having fl 0 = 0.3 or f l a = 0.1 does not yield very significant 
differences in the values of the long-string correlation length or velocity, it does produce very significant differences 
in the ratio of the long-string and loop densities to that of the background. 

Moreover, one should note that despite the significant drop in the number of loops produced (see||(a)), the ratio of 
the energy densities in loops and long strings decreases rather more slowly. This is because loops are slightly larger 
at formation and, since the average long-string velocity is decreasing, so is that of large enough loops (note that we 
assume that the RMS loop velocity at formation is equal to the RMS long stri ng v elocity at that time). Consequently, 
loops will live longer, since the redshift and gravitational radiation terms in (2.4) are velocity-dependent. 
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IV. CONCLUSIONS 


In this paper we presented the first discussion of cosmic string evolution in open universes, in the context of the 
generalized ‘one-scale’ model of Martins & Shellard |5|. We have shown that there is no linear scaling regime in an 
open universe, and that although the string density always decreases with respect to the critical density, it has been 
increasing relative to that of the background from z ~ , and it will become the main component of the universe 

sometime in the future. 

These differences with respect to the standard (flat universe) case only become significant fairly late in the matter- 
dominated epoch, so with respect to the string-seeded structure formation scenario we should only expect changes on 
very large scales—that one can easily estimate to be larger than the scales of the largest existing surveys. 

On the other hand, such large scales are of course relevant when one is comparing the cosmic microwave background 
anisotropies produced by cosmic strings to COBE data—thereby normalizing the string mass per unit length 11 
since this essentially involves an integration from the present time to the surface of last scattering. Thus the changes 
in the string network properties discussed in the present paper can significantly alter this normalization. This issue 
will be discussed in a forthcoming publication. 
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FIG. 1. Log-log plot of the scale factor a (relative to a(t eq ) = 1) as a function of cosmic time t (relative to t eq ) for a flat 
universe (solid line) and open universes with a present density S7 0 = 0.3 (dashed line) and 0, o = 0.1 (dotted line); both of these 
have h = 0.6 . 


5 



(a) 


1.6 

1.4 

1.2 

/ / 

/ / 

/ / 

CM 

1 

1 

2 / / 6 8 

CO 

o 

/ / 


/ / 

0.6 

- 

0.4 

/ 


(b) 


0.7 

-*x 


\\ 

0.6 

\r- 


\\ 

0.5 

\\ 


\ \ 

0.4 



\ \ 

0.3 

\ \ 


\ \ 


X 

-4 -2 

2 4 6'''-'.X 



2 


<d) 


0 


-2 


-4 


-4 



FIG. 2. Properties of a GUT long-string network in a flat universe (solid) and in open universes with a present density 
flo = 0.3 (dashed) and fi 0 = 0.1 (dotted), with h = 0.6 in both cases. Plots represent the ratio L/t (a), the RMS string velocity 
(b), and the log of the ratio of the long-string density to the critical (c) and the background (d) densities. The horizontal axis 
is labelled in terms of the logarithm of the scale factor (with a(t eq ) = 1); all plotted curves span the period between 10~ 10 t e9 
and 10 lo t e9 . 
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FIG. 3. Properties of the loop population of a GUT string network in a flat universe (solid) and in open universes with a 
present density = 0.3 (dashed) and fl 0 = 0.1 (dotted), both having h = 0.6. Plot (a) depicts the log of the number of loops 
produced per Hubble volume per Hubble time, while (b) shows the log of the ratio of the long-string (upper curve of each pair) 
and loop (lower curve of each pair) densities to the critical (b) and the background (c) densities. The horizontal axis is labelled 
in terms of the logarithm of the scale factor (with a(t eq ) = 1); all plotted curves span the period between 10^ 10 t eq and I0 10 t eq - 
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